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Abstract 


The  motion  of  elastic  firmed  projectiles  has  been  analyzed  by  various 
approximate  theories.  In  this  report  the  exact  equations  of  small  amplitude 
motion  are  derived  for  a  symmetric  missile.  The  aerodynamic  and  elastic 
syrrunetries  are  used  to  allow  the  use  of  complex  variables  to  describe  the  lateral 
motion  in  a  non-rotating  coordinate  system.  Although  the  resulting  equations 
are  both  ordinary  and  partial  differential  equations,  frequencies  and  damping 
rates  of  free  oscillations  are  obtained  from  an  ordinary  differential  equation  with 
boimdary  conditions.  Equations  for  a  permanently  deformed  bent  missile  are 
derived,  and  an  ordinary  differential  equation  for  the  forced  motion  of  a  bent 
missile  is  obtained.  Sample  calculations  for  a  finned  projectile  with  a  fineness 
ratio  of  20  show  resonant  motion  at  the  aerodynamic  frequency  as  well  as  at  each 
elastic  frequency.  The  nonlinear  roll  moment  associated  with  a  bent  missile  is 
computed  and  the  location  of  possible  spin-yaw  lock-in  is  determined.  The  flight 
motion  of  an  elastic  missile  is  shown  to  be  the  sum  of  two  elliptical  motions:  a 
low  frequency  pitching  motion  and  a  higher  frequency  flexing  motion.  The 
induced  drag  coefficients  for  both  motions  are  computed  as  functions  of  the 
missile's  elasticity. 
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1.  Introduction 


Modem  antitank  kinetic  energy  projectiles  have  become  very  long  finned 
missiles  carrying  very  dense  metal  rods.  Because  these  projectiles  have  to  resist 
very  large  aerodynamic  loads  due  to  their  high-velocity  flight  at  sea  level, 
designers  have  been  concerned  about  the  possibility  of  aeroelastic  deformations 
[1-7].  There  is  some  evidence  that  a  small  ntunber  of  these  projectiles  have  been 
forced  to  spin  at  rates  close  to  their  lowest  elastic  frequency  and  have  been 
subject  to  large  inelastic  deformations.  Mikhail  [3,  4]  and  Murphy  and 
Mermagen  [6, 7]  have  developed  special  solutions  that  showed  spin  lock-in  at  the 
lowest  elastic  frequency. 

Mikhail  [3,  4]  used  an  incorrect  expression  for  the  angular  momentum  and 
assumed  erroneously  that  a  spinning  eleistic  missile  could  perform  planar-flexing 
motion  and  a  planar-pitching  motion  in  a  coordinate  system  rotating  with  die 
missile.  His  numerical  calculations,  however,  showed  examples  of  spin  lock-in 
when  fin  damage  produces  a  roll-inducing  moment  sufficient  to  cause  a 
steady-state  spin  greater  than  the  lowest  elastic  frequency,  and  initial  spin  was 
zero.  Heddadj  et  al.  [5]  approximated  an  elastic  missile  by  two  rigid  bodies 
elastically  connected  and  studied  the  transient  motion  of  these  bodies  but  did  not 
consider  spin  lock-in. 

Murphy  and  Mermagen  [6,  7]  approximated  the  continuously  elastic  missile  by 
three  rigid  bodies  cormected  by  two  massless  elastic  beams  and  derived 
differential  equations  of  motion  for  this  three-body  system.  The  equations 
showed  that  it  is  impossible  to  cause  spin  lock-in  by  roll  inducing  moment  and 
zero  initial  spin  alone.  Murphy  [8,  9]  showed  that  the  combination  of  an  offset 
center  of  mass  (cm)  and  a  trim  pitch  moment  creates  a  nonlinear  roll  moment 
that  can  produce  spin-lock-in  at  the  aerodyricunic  frequency.  The  rear  beam  of 
the  3-body  theoretical  model  was  given  am  inelastic  lateral  deformation  and  spin 
lock-in  at  the  lowest  elastic  frequency  was  observed. 

Because  Mikhail's  [3,  4]  approximate  equations  are  far  from  being  theoretically 
correct,  their  results  can  be  dismissed.  Murphy  and  Mermagen's  [6, 7]  equations 
are  valid,  but  the  use  of  the  3-body  model  is  a  major  simplification  of  the  actual 
physical  problem.  It  is  the  purpose  of  this  report  to  derive  the  correct  equations 
for  a  continuously  elastic  symmetric  projectile.  This  derivation  exploits  the 
missile's  aerodynamic  and  elastic  symmetry  by  describing  the  transverse  motion 
with  complex  variables  and  using  a  nonrotating  coordinate  system  with  its  much 
simpler  equations  of  motion. 

Although  the  resulting  equations  form  a  combination  of  ordinary  and  partial 
differential  equations  containing  integrals  of  the  elastic  deformation,  free 
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oscillations  and  constant-spin  trim  motion  of  a  bent  missile  can  be  computed 
from  fourth-order  ordinary  differential  equations.  Induced  drag  for  various  free 
oscillations  are  deduced  as  quadratic  functions  of  their  amplitude.  Quadratic  roll 
moments  associated  with  a  bent  missile  are  formulated  for  trim  motion,  and 
equilibrium  spins  that  could  cause  spin-pitch  lock-in  are  located.  Niunerical 
integration  of  the  complete  set  of  partial  and  ordinary  differential  equations  for 
varying  spin  wovdd  be  necessary  to  verify  the  actual  occurrence  of  lock-in. 
Murphy  and  Mermagen's  3-body  theory  [6,  7]  would,  however,  be  a  good 
indicator  of  this  lock-in  behavior. 


2.  The  Coordinate  System 


If  we  consider  a  symmetric  missile  with  no  elastic  flexing,  we  can  define  a  non¬ 
spinning  coordinate  system  with  origin  at  the  rigid  missile's  cm  and  X-axis 
directed  eilong  the  axis  of  symmetry.  The  Z-axis  is  perpendicular  to  the  X-axis 
and  initially  downward  pointing,  whereas  the  Y-axis  is  determined  by  the  right 
hand  rule. 

The  elastic  missile  is  assumed  to  consist  of  a  very  heavy  elastic  circular  rod  of 
diameter,  d,  and  fineness  ratio,  L,  embedded  in  a  very  light  symmetric 
aerodynamic  outer  structure  that  may  be  longer  than  the  rod.  We  will 
conceptually  slice  the  missile  into  a  large  number  of  ttiin  disks  perpendicular  to 
the  X-axis  with  thickness,  dx.  For  the  elastic  rod,  an  elastic  nonspinning 
coordinate  system  is  defined  so  that  its  origin  at  the  center  of  the  central  disk  and 
the  X-axis  is  tangent  to  the  rod  at  this  point.  When  the  rod  flexes,  the  disks  shift 
laterally  perpendicular  to  the  X-axis  and  cant  to  be  perpendicular  to  the 
centerline  of  the  disks. 

The  angular  velocity  of  the  central  disk  in  these  nonspinning  coordinates  is 
defined  to  be  the  vector, 

a  =  {p,q,r).  (1) 

Because  the  coordinate  system  pitches  and  yaws  with  the  missile  but  does  not 
spin,  its  angular  velocity  vector  is 

£2  =  (0,q,r).  (2) 

A  dimensionless  vector  specifies  the  location  of  the  displaced  center  of  a  disk  at  X 
with  respect  to  the  origin  of  the  elastic  coordinate  system. 

8,(X,  t)  =  (X,Y,Z)d-'  =(x,5Hy,8,,).  (3) 

Because  of  die  mass  and  aerodynamic  symmetry,  it  is  convenient  to  represent 
transverse  displacements,  transverse  velocities,  transverse  angles,  and  transverse 
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angular  velocities  by  complex  quantities.  The  transverse  disk  displacement  and 
the  transverse  angular  velocity,  for  example,  are 

~  ^Ey 

Q  =  q  +  ir. 

The  lateral  location  of  the  cm  of  the  flexing  elastic  missile  can  be  expressed  in 
terms  of  die  missile's  local  mass  density  (x,  y,  z) , 


5c  (t)  =  (d/  m)|5E  (x,  t)>,  (x)dx ,  (5) 

where 

p,  (X)  =  JJp„ (x, y, z)dy(iz ,  and 

m  =  djp,(x)dx. 

Finally,  the  transverse  location  of  the  center  of  a  disk  relative  to  the  elastic 
projectile's  cm  is 


5(x,  t)  =  8E-6,, 


=  5y  +15^. 


(6) 


The  location  of  the  center  of  a  disk  relative  to  die  cm  and  its  velocity  relative  to 
die  cm  are  specified  by  the  vectors. 


5dc  =(x,5y,6j. 


(7) 

(8) 


The  transverse  components  of  this  relative  velocity  vector  cem  be  written  in  the 
form 


(9) 


where 


3.  Center  of  Mass  Motion 


The  velocity  of  a  disk  is 


(10) 
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where 


=  (v^ )  is  the  velcx:ity  of  the  cm. 


The  angle  of  attack,  a,  and  the  angle  of  sideslip,  p,  are  defined  at  the  central  disk 
for  the  elastic  missile.  If  the  transverse  velocity  components  for  this  disk  are 


=V4, 


(11) 


where 

^  =  P  +  ia,  and 


=  V, 


The  transverse  velocity  for  cm  of  the  elastic  missile  is 

Vy  +  iv^  =V4  +  S^d . 


(12) 


The  sum  of  equations  (9)  and  (12)  is,  therefore,  the  transverse  velocity  of  a  disk. 

+  M  +  ,,,, 

The  aerod)mamic  force  acting  on  a  rigid  projectile  is  assumed  to  be  proportional 
to  ^  and  Q  and  their  derivatives.  For  an  elastic  missile,  a  local  aerod)mamic 
force  needs  to  be  defined  in  a  similzir  manner.  The  x  coordinate  of  a  point  on  the 
rod  axis  lies  between  x,  =  — Z,/2  and  X2  =Z,/2.  The  aerod5mamic  structure 
incasing  the  rod  can  have  a  nose  windshield  of  length  Xjj  d  and  its  fins  could 
extend  beyond  the  end  of  the  rod  by  the  distance  Xq,  d  .  Thus  the  aerodynamic 
force  is  exerted  from  Xq  =  x,  —  Xm  to  Xj  =  X2  +  Xjj .  If  we  assume  the 
aerodynamic  structure  to  be  quite  rigid  beyond  the  ends  of  the  rod,  then, 

6E(x,t)  =  5E(x,,t)  +  (x-x,)r,(x„t)  ;  Xo<x<x,;  (14) 

8E(x,t)=6E(x2,t)+(x-X2)r2(x2,t)  ;  X2<X<X3.  (15) 


The  transverse  aerodynamic  force  exerted  on  a  disk  located  at  x  will  be  assumed 
to  have  the  form  appropriate  to  a  pointed  slender  body  [10, 11].  Thus,  the  force 
will  be  proportion^  to  the  local  angle  of  attack  and  its  modified  time  derivative. 
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where 


F  =  Fy+iF^,  and 
g,=pSV’l2. 

Cj-„  (oc)  are  two  normal  force  distribution  coefficients  per  axial  length  and  rj  is  the 

local  complex  angle  of  attack.  The  distribution  functions  derived  for  Mtmk's 
airship  theory  [10-12]  are,  for  example, 

dA 

=-2^5-';  Cf,=2AS-\  (17) 

where  A(a:)  is  the  local  cross-sectional  area. 

A  general  pointed  slender  body  relationship  between  Cyi  and  is 


JC 


(18) 


In  reference  [13],  computer  programs  for  calculating  c^j  for  finned  missiles  are 
given. 

The  complex  angle  between  the  missile's  centerline  described  by  S{x,t)  and  tine 
X-axis  win  be  denoted  as  F  and  is  the  spatial  derivative  of  d, 


as  ^  85^ 

dx  dx 


(19) 


The  local  angle  of  attack  is  the  angle  between  the  velocity  vector  and  the  surface 
of  the  missile,  which  is  assumed  to  be  parallel  to  the  centerline.  Because  the 
complex  angle  between  the  velocity  vector  and  the  X-axis  is  given  by 
equation  (13),  tiie  following  equation  results: 

7  =  ^  +  (4-ixe)b?/F-r.  (20) 

Equation  (20)  can  be  differentiated  and  the  small  angular  acceleration  terms  in  S 
and  Q  neglected  to  yield  a  relation  for  fi 

(21) 

The  total  transverse  aerodynamic  force  acting  on  the  missile  can  be  obtained  by 
integrating  equation  (16)  from  x^  to  Xj . 

F  =  IC„.^  +  C^A^Iv)+  iC^(Qdl¥)-J,  (()-  J,  (iXdMl  (22) 

where  Cfjj  are  listed  in  Appendix  A,  Jj{t)  are  listed  in  Appendix  B. 
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The  axial  aerodynamic  force  is  directed  along  the  X-axis  and  can  be 
approximated  by  use  of  an  axial  drag  coefficient  per  length  and  a  base-pressure 
drag  coefficient  , 

~  ~8i^dx  =  (^) 

where 

*3 

Cdx  =  fcD(x)ix  +  CDbp. 

Xo 

The  velocity  vector  of  the  cm  can  now  be  differentiated  as 

^  =  (n.v,.v.)+n*K.  (24) 

The  differential  equations  for  the  missile's  cm  are  therefore 

)  =  w  F  =  -g,Co ;  (25) 

/n(F|  +  F4  +  S^d-iVQ  )=F.  (26) 

Equations  (22),  (25),  and  (26)  can  be  combined  to  yield  a  simple  relation  between 
Q  and  ^ .  Small  terms  involving  and  €^,6  have  been  neglected. 

{V/dii-iQ)=-g^C,J  +  N,  (27) 

where 

C  =  r  -C 

^La  ^Not  ^D' 

S2  =gJnid,aDA 
N  is  defined  in  Appendix  B. 


4.  Angular  Motion  of  the  Projectile 


The  transverse  aerodynamic  moment  exerted  on  the  missile  is  computed  with 
respect  to  the  cm  and  has  contributions  from  both  the  transverse  aerodynamic 
force  and  the  axial  aerodynamic  force. 

dM  =-i(g^d%Cf{n-^Cf^{rjdlV))x-Co5]^,  (28) 

where 

M  =  My+iM^. 
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The  total  transverse  aerod5niamic  moment  acting  on  the  projectile  can  be 
obtained  by  integrating  equation  (28)  and  adding  ihe  base-pressure  moment. 

M^-iis,d1f„A*c^{idlv)*iC^iSdlv)-jM-jMdlv)-jM'  (2f) 

where  the  Cj^j  are  given  in  Appendix  A. 

Each  disk  pitches  and  yaws  with  angular  velocity  Q  +  zT  and  spins  with  spin 
rate  p^ .  The  spin  angular  momentum  about  the  center  of  a  disk  is  (d^/sjdm 
and  the  transverse  angular  momentum  about  the  center  of  a  disk  is 
{q  +  zf  jl^m .  Now  the  angular  momentum  of  a  disk  about  the  cm  of  the 
projectile  is  expressed  in  the  elastic  nonspinning  coordinate  system  and  the 
result  integrated  over  all  the  disks  to  yield  the  total  angular  momentum  of  the 
missile,  ht  d3mamics  terms,  we  will  neglect  the  mass  of  the  aerodynamic 
structure  and  assume  to  have  die  constant  value  of  miLd  .  Because  spin-yaw 
lock-in  is  caused  by  quadratic  terms,  quadratic  terms  are  retained  in  Pj  and  hj,^ , 


ii  =  (md 


(30) 


where 

p,  =p(i-rr/2)+R{Qf), 

=[2p(l-rT)+R{(Q-if^)l/l6,aiid 

iift+iifc=(Q+if+2pryi6. 

Equations  (7)  and  (8)  can  be  employed  to  yield  the  three  components  of  the 
angular  momentum  vector. 

H^=I^p  +  md^R\j,-Qj,};  (31) 

H,+iH,=l,Q-himd^(j,+jJ,  (32) 

where  Jjis  defined  in  Appendix  B. 

If  the  missile  is  assumed  to  fly  at  constant  pitch  angle  and  its  deformed  shape  is 
rotating  as  a  rigid  body  (Q  =  0,  ^  =  ip5 ),  the  axial  angular  momentum  is 

H,=hp,  (33) 


where 


i.=I.+(mdVL)][6|’-|rl7l6]dx. 

Mikhail  [3,  4]  computed  |<^(^)|^  for  a  flexing  missile,  neglected  the  |r|^  term  and 

assumed  that  equation  (33)  was  valid  for  the  axial  angular  momentum  for  any 
motion.  This  assumption  is  clearly  incompatible  with  equation  (31). 

The  derivative  of  the  angular  momentiun  vector  can  be  computed  in  the  usual 
way. 

^  =  [h,,H^,H,)+QxH.  (34) 

at 

For  a  slender  missUe  smd^L^/l2  and  jg  can  be  neglected  in  comparison  with 
jfg .  The  differential  equations  for  the  missile's  angtilar  motion  are 

IJ  +  (35) 

I,Q  -  ipI,Q  +  imd"  Jg  =  M .  (36) 

Equations  (27),  (29),  and  (36)  can  now  be  combined  to  give  a  simple  second  order 
differential  equation^  and  various  integrals  of  S . 

1,1  +  (a^  - ipl, )|  +  (a,  +  ipag )^  =  Je  +  Jn  +  md^  Jg,  (37) 

where 

^1  ““(8l<f)[CMo  ■•‘SsCLo^Mq  J' 

gj  =pSd/2m, 

g,  =I,/md^  ,and 

Je,  Jn  is  defined  in  Appendix  B. 

For  a  typical  missile  gj  is  less  than  10*^  and  the  second  term  in  a,  can  be 

neglected*  as  well  as  a  similar  term  in  Jn.  As  a  general  rule,  iQ  can  replace  |  in 
all  the  aerod5mamic  force  and  aerodynamic  moment  expressions.  The  rigid 
projectile  frequencies  and  damping  rates  can  be  computed  from  equation  (37)  for 


*For  lighter  than  air  dirigibles,  is  greater  than  0.1,  and  this  term  makes  a  statically  unstable 
airship  stable  in  flight  [14]. 


8 


the  right  hand  side  set  equal  to  zero.  According  to  Murphy  [15],  very  good 
approximations  for  the  frequencies  are 

K  =  k.  ±  ;  k  - 1,  2, 

=  +  {\±  pIjAI,(0^)pIj2I, 


(38) 


where 


=  is  the  rigid  projectile  zero-spin  frequency. 

Exact  relations  for  the  damping  rates  are 

^kR  ~  P^3  ~  x  )• 


(39) 


The  aerodynamic  roll  moment  is  the  X-component  of  the  aerodynamic  moment 
about  the  projectile's  cm.  The  linear  roll  moment  coefficient  for  a  rigid  finned 
projectile  is  usually  expressed  in  terms  of  a  roll-damping  coefficient  and  a  steady 
state  spin. 

The  steady-state  spin  is  usually  determined  by  a  differential  canting  of  the  fins 
caused  either  intentionally  by  the  designer  or  imintentionally  by  damage  to  the 
projectile. 

The  roll  moment  of  one  of  the  projectile's  thin  disks  is  the  sum  of  the  linear  roll 
moment  it  has  as  part  of  a  projectile  and  the  quadratic  roll  moment  induced  by 
the  transverse  aerod5mamic  force  acting  on  its  lateral  displacement  relative  to  the 
missile's  cm.  If  we  retain  only  the  dominant  Cj,  term  in  equation  (16),  the 
quadratic  roll  moment  has  the  form; 

=  (gi<*)criRKfe-r)+(8B  -ixQ)(d/V)^)dx. 

Therefore,  the  total  aerodynamic  roll  moment  acting  on  the  projectile  is 

M.=(g,dIC,),^+R{j,(t))],  (42) 

where  J,  is  defined  in  Appendix  B. 

The  nonlinear  roll  moment  from  equation  (42)  can  be  placed  in  the  spin  equation, 
equation  (35)  to  yield 

“Q  ^6  "82^9}“  82(^/)liiieiir‘ 

This  equation  is  nonlinear  due  to  the  retention  of  all  three  quadratic  terms. 
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5.  Flexing  Motion  of  the  Projectile 


Classical  beam  theory  for  a  circular  bezim  assumes  the  beam  to  be  slender  and 
the  elastic  moment  exerted  on  the  right  of  a  cross-sectional  disk,  M^,  is 

proportional  to  the  second  derivative  of  the  displacement  of  the  beam.  Although 
this  is  usually  stated  in  a  coordinate  system  rotating  with  the  beam,  the 
rotational  symmetry  of  the  beam  allows  us  to  state  the  proportionality  in  the 
non-spinning  coordinates. 

-iM.=(EI/d)^,  (44) 

OK 

where 

Me  +iM^r 

Id~*  =  JJ y^dydz  =  \\z^dydz  is  the  area  moment  of  rod,  and 
E  is  Young's  modulus  of  elasticity. 

A  similar  proportionality  for  the  elastic  shear  force  exerted  on  the  right  can  be 
obtained  from  beam  theory. 


F  =  — 


4  (El/dO- 


For  a  homogenous  circular  rod  with  constant  diameter,  El  is  constant  and  the 
elastic  shear  force  on  a  disk  is 

-^dx  =  -p,d  ©0  -T-f  ^ 

OK  OK 

where 

©^=El/p,d\ 

<&ohas  the  dimensions  of  a  frequency  and  appears  in  the  expression  for  the 
flexing  frequencies  of  a  free-free  beam  with  no  aerod)mamic  force  present. 

=(^a:/^)'®o'  (47) 

where 


=  4.730, 7.853,  10.996,  14.137... 
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A  very  common  beam-damping  assumption  is  the  Kelvin-Voight  model  [15], 
which  requires  the  beam-damping  shear  force  to  be  proportional  to  the  time 
derivative  of  tiie  elastic  shear  force.  This  time  derivative,  however,  must  be 
taken  in  rod-fixed  coordinates  that  rotate  with  the  rod. 

^dx  =  -p,d^a)lkc^,dx,  (48) 

dx 


where 


^dr  “ 


-1 


'd 

dt 

JJ 

''-I  ,  and 


A 

A:  is  a  small  dimensionless  beam-damping  coefficient. 

c*.  is  scaled  by  2o,  so  that  ^  =  1  corresponds  to  critical  damping  of  the  first 
elastic  mode. 


The  aerodynamic  shear  force  is  a  combination  of  the  aerodynamic  force 
distributions  Cf,(x),Cf2(x)  given  by  equation  (16)  and  the  drag-induced  moment 
distribution  Cd(x)5(x) appearing  in  equation  (28).  According  to  reference  [16], 
this  aerodynamic  shear  force  is 


dx 


dx  = 


dF 

dx 


5(c„S)- 


dx 


dx. 


(49) 


Although Cf,  and Cf 2  can  have  a  finite  number  of  jump  discontinuities,  Cp  must 
be  continuous  and  have  only  a  finite  munber  of  jump  discontinuities  in  its 
derivative. 

The  acceleration  of  a  disk  relative  to  the  cm  can  be  obtained  by  differentiating 
equation  (9).  The  linear  expression  for  the  lateral  components  of  this  acceleration 
is 

a  day  +  i<^dcz  =  •  (50) 

According  to  equation  (26),  the  lateral  component  of  the  acceleration  of  the  cm  is 
F/m ,  and  the  equation  of  motion  for  each  disk  is,  therefore, 

(p,cl')(8-ixQ  +  F/md)=^(F.  +F,  +F.).  pi) 
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dF 

The  aerodynamic  force  distribution  and  its  integral,  — ,F ,  are  available  from 

cbc 

equations  (16)  and  (22)  and  the  other  quantities  are  provided  by  equations  (46), 
(48),  and  (49). 


+  tOok  Ca4  -  gzL 


:„r  +  (c„f-c„8,)(d/V)-^^ 


E,^  +  E,4  +  ixQ-N 


(52) 


where 

E  j  (x)  is  defined  in  Appendix  B. 

Equation  (52)  is  based  on  the  usual  assumptions  of  neglecting  the  canting  inertia 
of  the  spuming  disks  and  any  shear  deformation.  These  small  effects  can  be 
calculated  if  desired  [17]. 

The  botmdary  conditions  at  x,  are  determined  by  the  aerodjmamic  force 
and  moment  exerted  on  the  overhanging  fins.  The  elastic  force  conditions  have 
an  additional  term  from  the  drag-induced  moment  distribution  [17]. 

^  V^’^^  +  kCa3(xi.t)+g4CD(xi)5(xi,t)=-(g4/g,)  dx  =  -g4f,;  (53) 

uX.  *0^^ 

^  +  kCd2 (xi , t)  =  (g4 /g, )  J |^(x  - X , - 6g,CD  j dx  =  g^m, ,  (54) 

where 

§4  =g2L/<Oo  /  and 

f),  mj  are  defined  in  Appendix  B. 

Similarly,  the  boxmdary  conditions  at  jCj  are  set  by  the  aerodynamic  force  and 
moment  exerted  on  the  nose  extension 


5x 

a^5E(x2,t) 


^  +  kc (x 2 , t) + g4CD (x2 )5(x 2 , t)  =  (g4 /g, )  J^dx  =  g4f2; 


(55) 


X  o  p* 

+  kCd2  (x 2 ,  t)  =  -(§4  /g, )  J  (x  -  x 2 )—  -  8g,Co 


,5x 

dx  =  -g4m2 .  (56) 


If  beam  damping  and  drag  eire  neglected  and  the  aerodynamic  structure  does  not 
extend  beyond  the  rod  (k  =  =  x^  =  X23  =  0 ),  the  boundary  conditions  reduce 
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to  those  for  a  free-free  beam. 


dx^ 


ac' 


=  0  ;x  =  Xi,jC2. 


(57) 


A  fourtih-order  partial  differential  equation  usually  has  fom-  boundary 
conditions.  In  this  case,  however,  there  are  two  more  conditions  at  x  =  0. 


5,(0)  =  0; 

a5,(o) 


ax 


=  0. 


(58) 

(59) 


These  conditions  look  like  die  usual  cantilever  boimdary  conditions  but  actually 
are  specified  by  the  motion  of  the  coordinate  system  which  is  attached  to  the 
missile  at  the  midpoint  and  tangent  to  it  there.  These  midpoint  conditions  can  be 
easily  satisfied  for  a  finite  element  calculation  and  need  not  be  used  as  boundary 
conditions.  Calculations  of  special  solutions  that  assume  a  single  frequency 
harmonic  time  vciriation  and  involve  ordinary  differential  equations  have 
difficulty  satisfying  the  midpoint  conditions.  In  this  report,  these  solutions  will 
be  calculated  by  pairs  of  differential  equations  in  x  for  the  aft  and  fore  section  of 
the  rod  (xi  <  x  <  0;  0<  x  <  X2).  Equations  (53),  (54),  (58),  and  (59)  are  boimdary 
conditions  for  tiie  aft  part  and  equations  (55),  (56),  (58),  and  (59)  are  boimdary 
conditions  for  the  fore  part.  This  use  of  six  boundary  conditions  implies 

discontinuities  in - f-, - at  the  central  junction  point.  These  six  boundary 

dx  3x 


conditions  for  the  ordinary  differential  equations  in  x  are  not,  however, 
independent,  and  it  can  be  shown  that  equations  (53)  and  (54)  imply  equations 
(55)  and  (56)  and  vice  versa  (Appendix  C).  Thus,  the  calculations  should  not 
show  discontinuities  in  shear  force  and  moment  at  the  junction  point. 


6.  Bent  Missile 


If  the  rod  were  inelastic2dly  deformed  on  launch,  its  shape  would  be  represented 
by  the  sum  of  a  fixed  deformation  rotating  with  the  missile  and  an  elastic 
deformation, 

where 

8EB(0)=%^=0.and 

dx 
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The  local  inclination  of  the  rod  and  the  lateral  location  of  the  cm  become 


(61) 

(62) 


r=f+r,c"; 

5.=5.+5«e*, 

where 

~  dL 

OX 

■p  _ 

dx  ' 

6,  =(l/L)j6Edx,and 

^cB  ~  (V^)  JSeb^  . 

*1 

The  location  of  the  centerline  of  the  disks  relative  to  the  projectile  cm  is 

S  =  S,-l+S,e‘*, 

where 


For  the  rigid  aerod)mamic  structure  extension, 

^EB  (■*)  ~  ^EB  (■*'1  ~  (^1 )  '  Xq  ^  X  ^  X^  , 

^EB  W  “  ^EB  (^2  )■*■  ~  ^2  (-^2  )  '  X2  ^  X  <  X^  . 


(63) 


(64) 


Now  produces  2in  aerodynamic  force  and  moment  that  rotates  with  the 
projectile.  If  the  fins  are  bent  with  respect  to  the  rod,  an  additioncil  differential 
force  that  rotates  with  the  missile  is  produced  and  must  be  added  to 
equation  (16). 

^^BF  =giCfjrB,r(xy*dx,  (65) 


where 

I’bf  W  ®  measure  of  bent  fin  damage. 

The  aerod5mamic  force  acting  on  the  bent  missile  is 

F  =  -g,Lc».4  +  (CNd  +C^,)(5d/V)-J,(t)- j,(tXd/V)-C^,e“J,  (66) 

where  Cnbf  is  defined  in  Appendix  A. 
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If  the  fin  damage  extends  beyond  the  rear  end  of  the  rod,  it  will  exert  an 
additional  force  and  moment  on  the  rod  at  x  =  Xy.  This  bent  fin  force  term  also 
has  to  be  included  in  die  quadratic  roll  moment  appearing  in  equation  (43). 

8xP‘^^{j7~QJ6“S2(j9'*'  ^9BF  )}  “  §2  )linear  • 

where 

0  _ 

J9BF  ~  ^  JcfjrBpbdx . 

*0 

Equation  (27)  for  the  cm  of  the  elastic  missile  becomes 

(V/d)(4-iQ)=-ga(Cu.5-C*,e'*)+N.  (68) 

The  aerodynamic  moment  as  given  by  equation  (29)  becomes 

M  =  -ife,d)Lc^^4  +  (c^  +C„,X4d/v)-J,(t)- j.(tXd/V)-J,(t)-C„,,e‘‘J,  (69) 

where 

Cmbf  is  defined  in  Appendix  A. 

The  projectile  angular  motion  equation  (36)  becomes 

I.Q-ipI.Q— i(g,d)[c^4+(CK„  +C„,)(4d/v)-C„„e'*]-iJ,  -imd%.  po) 


In  file  flexing  motion  equation  (51),  S^iis  replaced  by  in  the  elastic  and 
damping  force  terms,  and  the  bent  fin  terms  are  added  to  F  and  Fa  to  obtain  the 
partial  difierential  equation  for  a  bent  missUe. 

+ ®o  -  g^Lr  c„r + |:(cf2r  - 

dt  dx  ot  ox 


(71) 


=  E,4  +  E24  +  ixQ- 


N-E^pe^ 


In  the  boundary  conditions  specified  by  equations  (53)-(56),  is  replaced  by 
Sg  in  the  partial  derivatives  on  the  left  side  of  the  equations.  The  conditions  at 
Xj  £ire  modified  by  force  or  moment  contributions  from  any  fin  damage  between 
Xq  and  X,  that  may  have  occurred. 


7.  Trim  Solution 


The  general  solution  for  the  motion  of  a  bent  missile  requires  the  numerical 
integration  of  one  partial  differential  equation  (equation  [71])  and  three  ordinary 
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differential  equations  (equations  [67-69]).  A  simple  special  solution  is  that  for 
the  steady  state  motion  of  a  bent  missile  with  a  constant  spin.  This  motion  is 
called  trim  motion  and  has  the  form 


II 

(72) 

II 

(73) 

(74) 

11 

+ 

II 

(75) 

r=(rB+fT^^* 

(76) 

For  this  motion,  beam  damping  is  zero  ( =  0  )• 

Equation  (68)  can  now  be  used  to  eliminate  Q  and  Q  from  equations  (70)  and  (71) 
and  equations  (72-76)  inserted  to  3deld: 

bj^T  “  Jet  JnT  JfiF  P  *^67  ' 

(77) 

-  E.  + E,5^  -  g  A.  ^  =  E,5,  -  N;  - E-„  +  E„ 

dx  dx  dx 

(78) 

where  JCT,JNT.JBF.b,,Ej,EB,N;,EBF  are  defined  in  Appendix  D. 

Boundary  conditions  as  given  by  modified  equations  (53)-(56),  (58),  amd  (59) 
become 

~et(  *^  +  g^cD(x,)6T(x,)  =  -g4f,T; 

dx 

(79) 

''  =  g4m,T; 

dx 

(80) 

"^  +  g4CD(x2)5T(x2)=g4f2T; 
dx 

(81) 

dx 

(82) 

6^(0)  =  0; 

(83) 

dd^ri^)  _  n 

dx 

(84) 

where  fjr.nijT  are  defined  in  Appendix  D. 
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These  boundary  conditions  can  be  satisfied  by  expressing  ^ET  as  the  stun  of 
three  functions  of  x, 

^£tW=WiW+^2>»'2W+^3>*'3W-  (85) 

The  first  function  is  ihe  solution  of  the  following  differential  equation  and  initial 
conditions 


^-E,^+E,w,  -N;  -E;,  +E, 


dx^ 


dx 


dx 


(86) 


where 


yv, 


fn)  =  _  ^^^^1  (o)  ^  (Q)  _  0 

^  dx  dx^  dx^ 


The  otiier  two  functions  are  solutions  of  the  homogenous  equation  with  different 
initial  conditions. 


dV„ 

dx" 


--E. 


dw„ 

dx 


■  +  E5W_=0;  m  =  2,3. 


(87) 


where 


,^(0)  =  ^  =  £^  =  0,and 
^  dx  dx' 

'w^Co)  d'w3(0) 

dx'  dx' 


Equations  (86)  and  (87)  can  be  integrated  from  zero  toward  negative  x  or  toward 
positive  X.  w„,(x)  denotes  for  negative  x,  and  w„2(x)  denotes  for 
positive  X.  are  values  of  which  satisfy  equations  (79)  and  (80), 

andfi22»-®32  those  values  which  satisfy  equations  (81)  and  (82).  The  values 

should  be  close  to  each  other  and  their  difference  is  a  measure  of  the  accuracy  of 
various  approximations  and  the  numerical  integration  process. 

The  presence  of  the  constants  in  equations  (77)  and  (86)  and  the 

bormdary  conditions  present  some  difficulty  because  they  involve  integrals  of 
and  its  first  derivative.  Our  numericeil  technique  is  to  replace  these 
constants  witii  guessed  values  Nxo,4to>5co  use  an  iterative  three- 

dimensional  Newton's  Method  to  converge  on  the  proper  values.  The  partial 
derivatives  used  in  Newton' s  Method  are  computed  by  ttie  Chapman-Kirk  [18] 
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technique.  If  drag  is  neglected,  the  need  to  know  6^  vanishes,  and  this  reduces 
to  a  two-dimensional  process.  This  process,  of  course,  depends  on  reasonably 
good  first  guesses. 

For  trim  motion  p  =  j7  =  0,Q  =  ip^^^  and  equation  (67)  become  a  simple 
equality  of  two  functions  of  p. 

fi(p)  =  fM.  (88) 

where 

fl=P-Pss> 

f,  =-R|(pV8a)lTj.T  +J,T)(c»d/v)-'  ,and 

JjT  Jcf|[(^,  — r,  — rpp)+i(8CT  —  )^>d/ V)]5,dx. 

Xo 

Equilibrium  values  of  spin  are  determined  by  the  intersection  of  these  two 
curves.  The  stability  of  an  equilibrium  spin  can  be  foimd  by  the  integration  of 
equations  (67),  (68),  (70),  and  (71),  for  nonconstant  spin  near  an  equilibrium  spin. 
The  stable  equilibrium  spins  locate  spin  lock-in  possibilities. 


8.  Transient  Solutions 


A  rigid  S5ntnmetric-finned  missile  flying  with  constant  spin  has  two  natural 
frequencies,  and  ^2r  >  ''vhere  =  -^2r  •  Each  of  the  flexing  frequencies 
given  by  equation  (47)  would  give  rise  to  two  coning  frequencies,  ±  .  The 

frequencies  present  in  the  motion  of  an  elcistic  projectile  would  form  an  infinite 
sequence  where  the  first  two  would  be  related  to  and  (f>2g ,  whereas  dte  later 
ones  would  evolve  from  i-COg,  i.e.,  {i>2K+\=^K  ■>^ik+2  =  ~^k) ' 

numbered  modes  have  positive  frequencies  and  are  called  positive  modes, 
whereas  the  even  numbered  modes  are  called  negative  modes. 

Transient  solutions  of  equations  (70)  and  (71)  are  solutions  of  the  homogeneous 
form  of  these  equations,  (dg  =  Fg  =  C^p  =  C^^gp  =  O).  Special  transient  solutions 
of  equations  (70)  and  (71)  have  the  form 

^  =  (89) 

8E=>('i(xKe‘-'.  (90) 

where. 
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A 

T  =  1C 

Equation  (68)  can  be  used  again  to  eliminate  Q  and  Q  from  equations  (70),  (71), 
(89),  and  (90)  inserted  to  peld 

^Ik  “  ^Ek  ■'■^Nk  -^k^Sk'  (^1) 

^  ^  ~ 

where  b ,  Jg^ ,  ,  Ejk  >  Vkc  >  Nr  are  defined  in  Appendix  E. 

These  equations  and  their  boundary  conditions  can  be  solved  in  a  manner  quite 
similar  to  that  used  for  the  trim  solution  by  the  introduction  of  three  auxiliary 
functions. 

W=  1^4  W+^5>V5  W+^6>V6W-  (93) 

W4  is  the  solution  of  the  following  differential  equation  and  initial  conditions: 

E»  ^  -  E„  E,^w.  =  -E„  -  n;  +  ^ .  (94) 

where 

dx  dx^  dx^ 

The  other  two  functions  satisfy  the  homogeneous  equation  with  different  initial 
conditions 


where 


"ek 


d‘’w  dw 

^  vv  m  p  m 

4  ^4k  j 

dx  dx 


+  E5kW„  =0;  m  =  5,6. 


w. 


w. 


dx  dx^ 
dx^ 


dx 

d^w^i^)  d^w^ 


=  0, 

=  0/and 


(95) 


The  parameters  are  values  of  B^,B^  that  satisfy  the  boundary  conditions 

at  X,: 
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(96) 


^7^  +  g4E;lco(\i/k  -ykc)  =  -g4E;[f,k; 

dx 

^^  =  g4E-Xk-  (97) 

dx 

Similarly,  B  52 ,653  perform  the  same  function  at  Xj . 

+  g4E;LcD  (Vk  -  M^kc )  =  g4E6lf2k ;  (98) 

dx 

^r^  =  -g4E;!cm2k.  (99) 

dx 

The  presence  of  the  parameters  complicates  these  differential 

equations.  involves  A^and  N|j.  The  definition  of  involves  integrals  of 
\|;^  and  its  first  derivative  and  must  satisfy  equation  (91)  that  contains 
and  integrals  V|/^  and  its  first  derivative.  Guessed  values  of  these  parameters 
Nko,Ako,V|/kcoare  used  in  equations  (94)-(99)  and  another  three-dimensional 
Newton's  Mettiod  is  used  as  well.  For  frequencies  greater  than  3©^ ,  a  good 
estimate  for  is  —  Nq Aj^J, . 


9.  Flight  Motion 


The  complete  flight  motion  of  an  elastic  bent  symmetric  projectile  with  constant 
spin  is  described  by  the  sum  of  a  trim  mode  and  an  infinity  of  transient  modes. 
If  we  limit  ourselves  to  tite  aerodynamic  modes  and  the  first  two  elastic  modes, 
the  pitching  motion  and  the  lateral  motion  of  the  forward  end  of  the  rod  can  be 
expressed  by  two  sums  of  five  terms. 

4  =  +  T,e'''‘  +  ;  (100) 

S,(x„t)  =  S^(x,)e*'  +t,e*''  +t,e*’'  +t,e*’'  +t.e*'',  (101) 


where 

Tk  =¥k(x2)Tk  ='l/k2Tk- 

Because  \|/,2  and  \\I22  usually  less  than  0.5  and  the  cingular  motion  is  usually 
less  than  0.1  radian,  the  motion  described  by  the  aerodynamic  modes  is  primarily 
a  pitching  motion  and  is  best  characterized  by  (K,,K2)l  |\|;32|and|\j/42|  are. 
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however,  usually  greater  than  fifty  and  the  motion  described  by  these  flexing 
modes  is  primarily  a  flexing  motion.  For  this  reason  we  will  characterize  this 

flexing  motion  by  Ihe  modal  amplitudes  of  tfie  rod  lateral  motion  (k3o,K4o). 

For  a  nonspinning  missile,  equations  (91)  and  (92)  imply  an  important  pairing  of 
solutions.  If  is  a  solution,  A^.vj/kis  a  solution. 


Az  =  A,  M/j  =  v|/, 

A4=A3  X|/4=V3. 


(102) 


Thus  for  zero  spin,  the  negative  modes  have  the  same  dzunping  as  the  positive 
modes,  and  tixeir  waveforms  are  conjugates  of  the  positive  waveforms.  For  small 
spin,  the  magnitudes  of  the  waveforms  are  essentially  equal, 

kl2|  =  k22tk32|  =  k42l-  (103) 


The  transient  angular  motion  and  the  rod  lateral  motion  are  botii  expressed  by 
the  sum  of  two  damped  elliptical  motions  with  frequencies,  •  The  initial 

semi-major  axes  of  the  angular  modes  are  Kj^  +  aiKl  (K.„+K«,)|v,2r. 
whereas  the  initial  semi-minor  axes  of  the  angular  modes  are  |K,o  “Kjoj  and 
|\|/32|'' .  Similarly,  the  initial  semi-major  axes  of  the  rod  lateral'modes 
(Kjo  +K2o)|Vi2|  aiid  (^30  +K40),  whereas  the  initial  semi-minor  axes  are 


K30-K40 


are 
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and  IK30-K 


“■40 


If  either  motion  goes  througji  zero,  the  ellipticcil  motion  would  be  planar.  For 
gun  laimch,  both  the  angular  and  the  flexing  motion  are  forced  to  be  near  zero 
initially,  and  "near  planar  motion"  is  expected  after  launch.  Guidos  et  al.  [19] 
have  analyzed  the  flight  motion  of  long  elastic  finned  projectiles  in  the 
U.S.  Army  Research  Laboratory  (ARL)  Transonic  Experimental  Facility's  non- 
rotating  coordinates  and  have  observed  near  planar  nose-tip  motion  with  an 
amplitude  in  excess  of  0.5  diameter. 


10.  Induced  Drag 


The  aerod3mamic  drag  force  is  directed  along  the  velocity  vector  and  the  linear 
normal  force  on  a  disk  is  perpendicular  to  the  disk  axis,  which  is  canted  at  an 
angle  of  ^  “  T  with  respect  to  the  velocity  vector.  Thus,  the  normal  force  has  a 
quadratic  contribution  to  the  drag.  This  contribution  is  called  induced  drag  and 
can  be  computed  for  a  pitching  and  flexing  missile  as 
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(104) 


Sl^DI 


dx . 


Although  the  nonlinear  axial  flow  along  the  rod  has  a  quadratic  contribution  to 
drag,  this  contribution  is  usually  significantly  less  than  the  induced  drag. 

The  normal  force  distribution  will  be  approximated  by  its  dominant  term, 

Cnte-r). 


where 


CD,  =  Kfe-rX5-r)ix 

*0 


(105) 


)  - 

Jd  =  Jc„rrdx. 

Xo 

The  first  term  in  equation  (105)  is  the  weU-known  expression  for  the  induced 
drag  of  a  rigid  missile.  The  induced  drag  for  a  particular  angular  mode  is 

^DIk~^k^kO' 


where 

^k  ~  ^Na  —  ^Ik  “  ^Ik  ^Dk  ' 


'Dk 


=  JCf, 


dv|/k  dvj/^ 


*0 


dx  dx 


dx. 


The  average  induced  drag  for  a  missile  whose  motion  contains  aU  five  modes  of 
equations  (100)  and  (101)  can  be  computed  from  the  svun  of  the  pure-mode 
induced  drag  contributions  and  the  trim  induced  drag. 


C„,  =C,(Kf.  +K^.)+C,(k;.  +Ki,)+CB„, 

where 

C3  ~C3|\|/33|  ,  and 


(107) 


^^3 

^DIT  ~  “^t|  ^ 


Xo 
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11.  Numerical  Results 


In  references  [6]  and  [7],  numerical  results  were  obtained  for  two  specific 
missiles.  The  first  was  a  flare-stabilized  25-cal.  rod  flying  at  18,000  ft/s 
considered  by  Flatus  [1]  and  the  second  was  a  fin-stabilized  20-cal.  rod  flying  at 
6000  ft/s.  The  flare  on  the  Flatus  missile  was  7.5-cal.  long  and  the  flare  base 
diameter  was  4  cal.,  i.e.,  4  rod  diameters.  The  finned  missile  has  a  1-cal.  nose 
extension  and  a  1-cal.  fin  extension  (Figure  1).  The  necessary  parameters  for 
these  missiles  are  given  in  Appendix  F  and  converged  values  of  for 

the  finned  missile  are  given  in  Table  1.  The  slender  missile  values  of  the 
moments  of  inertia  are  given  in  parenthesis  are  used  in  the  nximerical  calculation 
of  the  transverse  motion.  The  actual  value  of  h  should  be  used  in  an  integration 
of  the  roll  equation  (67). 

A  measure  of  the  flight  flexibility  of  a  missile  is  the  ratio  of  the  first  elastic 
frequency  to  the  rigid  missile  aerodynamic  frequency,  (j  =  a)^l(Og.  The 
equations  for  the  transient  solution  derived  in  section  8  have  been  solved  for  the 
first  aerodynamic  frequency,  and  the  results  plotted  versus  <t  for  both 
missiles.  In  Figure  2,  it  is  shown  that  this  frequency  for  the  flare-stabilized 
missile  is  98%  of  the  rigid  missile  value  for  cr  =  10  and  decreases  as  <7  decreases. 
We  see  that  the  theory  of  this  report  falls  between  the  3-body  theory  and  the 
Flatus  theory. 

The  remainder  of  the  numerical  results  was  obtained  for  the  finned  missile.  In 
Figure  3,  tiiis  report's  theory  is  compared  with  the  3-body  theory,  and  we  see 
that  it  has  a  similar  variation  but  predicts  lower  frequencies.  Note  that  the 
frequency  for  cr  =  5  is  neeir  60%  the  rigid  missile  value.  The  first  two  positive 
flexing  frequencies  are  plotted  versus  ct  in  Figure  4.  For  small  o 

both  frequencies  are  slightly  larger  than  the  free-free  beam  values,  but  they 
approach  these  values  as  ct  grows. 

The  magnitude  of  the  forward  end  motion  is  plotted  against  o  in  Figure  5.  At 
cr  =  5,  we  see  that  0.1  radian  of  the  positive  aerodynamic  mode  causes 
0.18  diameter  of  flexing  motion.  For  the  flexing  frequencies,  one  diameter  of 
flexing  motion  is  associated  with  0.018  radian  (1.0°)  of  angular  motion  for  the 
first  positive  mode  and  0.072  radian  (4.1°)  of  angular  motion  for  the  second 
positive  mode. 

The  forward  end  of  the  rod  bends  away  from  the  velocity  vector  for  an 
aerodynamic  mode  (xi/j2  ^  O)  while  it  bends  toward  the  velocity  vector  for  both 
the  first  positive  symmetric  mode  and  the  first  positive  antisymmetric  mode. 
The  flexing  shape  of  the  entire  rod  is  very  close  to  planctr  for  each  mode.  Thus 


R{v|/^  (x)/\|/^2}  describes  the  shape  of  the  rod  during  a  modal  motion  and  is 
plotted  versus  x  for  k  =  1, 3, 5  in  Figure  6.  For  the  first  two  positive  frequencies, 
the  shape  is  quite  similar  and  the  rear  end  of  the  rod  is  deflected  about  15%  less 
than  the  forward  end.  The  third  positive  mode  is  antisymmetric  and  rear  end 
deflects  about  5%  more  than  the  front  end  and  in  tiie  opposite  direction. 

The  damping  rates  of  the  three  modes  for  no  beam  damping  (k  =  o)  are  plotted 
versus  a  in  Figure  7.  For  cr  =  5 ,  the  first  aerodynamic  damping  rate  is  85%  of  the 
rigid  body  damping  rate,  while  the  rates  for  the  second  and  third  positive  mode 
damping  rate  are  75%  emd  35%,  respectively.  In  Figure  8  X,j/cOr  is  plotted 

versus  p/co^  for  two  non-zero  values  of  beaun  damping  and  a  =  5 .  for 

this  value  of  a  is  0.6  and  beam  damping  reduces  the  damping  rate  when  p/to^ 
exceeds  this  value.  The  behavior  was  first  observed  Flatus  [1]  and  predicted  by 
the  3-body  theory  [6, 7]. 

The  induced  drag  coefficient  for  the  aerodynamic  mode,  Cj/C^, ,  is  plotted 

versus  a  in  Figure  9.  It  is  12%  less  than  that  for  a  rigid  middle  at  o  =  5. 
Although  the  nose  of  the  elastic  missile  bends  to  a  larger  angle  of  attack,  the  fins 
bend  to  a  lesser  angle  of  attack  and  the  net  effect  is  less  induced  drag.  The 
induced  drag  coefficients  for  the  two  flexing  modes  are  plotted  versus  a  in 
Figme  10.  Here  we  see  that  for  the  same  amplitude  motion  the  induced  drag  for 
the  symmetric  mode  is  more  than  50%  greater  fiian  that  for  the  antisymmetric 
mode.  Because  both  ends  of  the  missile  bend  to  a  lesser  angle  of  attack  for  the 
antisymmetric  mode,  this  mode  should  have  much  smaller  induced  drag. 

For  the  bent  missile  calculations  we  wiU  asstune  no  fin  damage  (^BF  =  O)  and  the 
bent  rod  wiU  be  described  by  a  pair  of  quartic  curves. 


8^  =d^^x^ +d2xX*  -10<x<0 
=  +  dj^x*  0  <  X  <  10 

Equation  (62)  provides  a  relation  for  5^8 . 

5.,  =(d„  +d„XL/2)76  +  (d3,  +d^)(L/2)VlO. 


(108) 


(109) 


We  wiU  assxune  the  rear  of  the  rod  is  imdeformed  d,,  =  r/,2  =  0  and  tiie  values  of 
dj2,d22  are  given  in  Appendix  F. 

In  Figtue  11,  the  magnitude  of  the  ratio  of  the  trim  angle  of  attack  to  the  zero 
spin  trim  angle  of  a  rigid  missile  is  plotted  versus  spin  for  c  =  5 .  The  zero  spin 
trim  angle  of  the  elastic  missile  is  about  four  times  that  of  a  rigid  missile.  The 
first  resonant  peak  is  near  p/oJr  =  0.6  cind  shows  a  five-fold  magnification  of  the 
zero  spin  elastic  trim.  A  smaUer  resoiumt  peak  occurs  at  the  elastic  frequency 
near  5  and  much  larger  resonant  pecdc  near  14.  The  magnitude  of  the  rod 
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deflection  is  plotted  versus  p/(0r  in  Figure  12.  The  trim  deflection  for  the  first 
elastic  mode  is  about  eight  times  that  for  the  aerod)mamic  mode  while  the  trim 
deflection  for  the  second  elastic  mode  is  three  times  that  for  the  first  elastic  mode. 

In  Figure  13,  fa  is  plotted  versus  p/cOr  .  A  sample  fi  line  for  Pss/wr  =  6.5  is  also 
shown.  The  intersections  of  these  curves  are  loci  of  possible  spin-yaw  lock-ins. 
The  actual  occurrence  of  lock-in  would  be  shown  by  integration  of  the  partial 
differential  equation  for  flexing  motion. 


12.  Summary 


1.  A  system  of  equations  has  been  derived  for  the  motion  of  a  symmetric  elastic 
missile.  This  system  consists  of  three  ordinary  differential  equations  and  one 
fourth-order  partial  differential  equation.  Each  of  these  equations  contains 
integrals  of  the  elastic  deformation  of  the  missile. 

2.  A  fourth-order  ordinary  equation  with  boundary  conditions  has  been 
obtained  to  determine  natural  frequencies  and  damping  rates  of  the  motion. 

3.  A  simple  special  ordinary  differential  equation  for  the  trim  motion  of  a  bent 
missile  has  been  obtained. 

4.  Equations  for  nonlinear  roll  moments  for  a  bent  missile  have  been  obtained 
and  a  condition  for  possible  spins-yaw  lock-in  stated. 

5.  Numerical  results  have  been  obtained  for  the  natural  frequencies,  flexing 
waveforms,  damping  rates,  and  induced  drag  coefficients. 
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Figure  2.  for  flare-stabilized  rod, 


Figures.  ^x/cOd  vs.  <T  for  filmed  missile. 
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Appendix  A.  Integrals 
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Appendix  6.  Functions 


Ji(t)=  /c„rdx 

Xo 

13(1)=  JCfjFxdx 

Xo 

J4(t)=  j(Cf2r-Cf,5E)x(iX 
xo 

J,(t)=iTc„[fe-r)+(8H  -ixQ)(d/V)]5dx 

J.(t)=-(g4j.+j4(<l/V)+J5] 

JN(t)=[l,N-(ipI.)N]('i/V) 


h(^)=  /cD5dx  +  8(x,)cDbp 
J6(t)=(VL)/x6Edx 
J,(t)=-(i/L)?^5  +  f„f)lx 
J,(t)=-(l/16L)?(f-2ipr)b 


c*  =  2(of 


a' 

ax^ 

-ipSi 


N(t)=g.[j,+j,(d/V)]-8. 

E,W=gj(C„.-LCn) 

E.W=g=(CN*  +C„,  -Lc„)(d/V) 

Ebf  “SlC^NBF  “ECfiPpp) 
c„(x)  =  2c,j-xc„ 

f„(x,t)  =  {f-2ipr  +  iQVl6 

f.(t)=I,fe(t)-r(x„t)]+|l„4(t)+I,8B(x„t)+(l,-I,X'(x„t)J(d/v) 

f,(t)=ij5(t)-r(x„t)]+[i„4(0+M,(x,.t)+(i4-i.)(-(x„t)](d/v) 

in,(t)=I,[5(t)-r(x„t)]+[l„4{t)+lA(x,.t)+(l.-l7Xtx..‘)kv)-ni,D 

m7(l)=I.fe(0-r(x7,t)]+[l„4(t)+IA(x„t)+(l„-I.)f(x„t)](d/V)-in, 

^ID  (0  ”  ^ID  [^E  (^1  >  0 ““  ^3D^(^1  >  0 

^2D (0  ~  ^2D  [^E  (^2  » (0]“^  ^4D^(^2  »  0 
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Appendix  C.  Relationship  Between  Boundary  Conditions 


If  we  make  the  slender  missile  assumptions  (I,  =  md^L^/l2 ;  I^/l,  sis  As  =  0), 
the  boundary  conditions  at  x,  and  are  directly  related.  This  can  be  shown  by 
first  integrating  equation  (51)  from  Xj  to  x  at  fixed  time. 

l8dx-iQ(x= -x;)/2-g,[c^^-J, +[(C„ +C^X- j!](d/V)](x->i,)/L  = 

*1 

,ra^5E(x)  a^8E(xi)l  >  ,r  /  x_  /  n 

”^0  0^3  0^3  ^®ol^d3W  ^d3V^l/J  (C-1) 

-g,  |[c„fe-r)+t:„{8.  -xO+Cr2(24-r)td/V)]dx 

-g2[cD(x)5(A-CD(xi)5(x,)] . 

Equation  (C-1)  is  now  evaluated  at  Xj  using  the  integrals  of  C.l. 

^  +  kc^ (x2)=  +  kc^j (x,) - g4 [cd (xj )5(x2 )- Cd (xi )5(x, )] 

+  84(^1  +^2)  ' 

where  f  j  is  defined  in  Appendix  B. 

According  to  equation  (C-2),  boundary  condition  (53)  implies  boimdary 
condition  (55). 

If  equation  (C-1)  is  integrated  from  Xj  to  Xj ,  five  single  integrals  and  six  double 
integrals  result.  Each  of  the  double  integrals  can  be  simplified  by  integration  by 
parts  and  are  given  in  section  C.3. 

Under  our  stated  assumptions,  equation  (36)  becomes 

(LVl2)Q  +  iJs  =-ig2[cMa^- J3  +  ((CMd  +Cm,X- Akv)- J5] .  (C-3) 

Equation  (C-1)  is  now  integrated  and  simplified  by  use  of  sections  C.1-C.4  and 
equations  (53)  and  (A-3), 

+  kca2 (x2 )  =  ^  +  kCi2 (x, )+  g4(mi  +  )  ,  (C-4) 

Sx  ox. 

where  mj  is  defined  in  Appendix  B. 

Boundary  condition  (54)  is  equivalent,  therefore,  to  boimdary  condition  (56). 


C.1  Integrals  in  Equation  (C-2) 


|sdx  =  0  Jcf,dx  =  CN„ -I, -I2 

l|  X, 

JcfjFdx  =  Jj  —  IjE,  —  IjFj  Jcfjdx  =  —  Ij  —  Ij 

‘1  *1 

}(Cf2  -xcn>ix  =CNq  -(l„  -l5)-(l,2  -Ifi) 

^1 

~  “  “'^2  ~  ^l^EI  —  (^3  ~  ^5  X'l  ~  I2^E2  ~  (^4  ~  ^6  K'2 


C.2  Integrals  in  Equation  (C-4) 

(1/2)  J(x^  -  xf  )dx  =  -  L^/12 


L”'  J(x-x,)dx  =  L/2 


^2 

\ 


0^6e(x)  a^5E(x,) 


_^'5e(x^)  a"5E(xi)  g^5E(xi). 

ax"  ax"  ax" 


ax"  ax" 

J[Cd3  (x)-  Cd3(x,  )]dx  =  Cd2(x2 )-  Cd2  (x, )-  Cd3  (x,  )L 

^2 

JcD8dx  =  Jj-l3D-I 


■^40 


C.3  Multiple  Integrals  in  Equation  (C-4) 


Xjx 


J  Jadxdx  = -LJg 

ii  %x 

|j[c„dxdx  =  L[(l/2)C„.  -I,]-[C„.  -I,  -I,] 

l,Xi 

I  J[c„rdx  dx  =  l[(i/2)i,  -  I.r.  ]-  [j,  -  ijF,  -  i,rj 

I,  X, 

I  Jc^dxdx  =  L[(I/2)C«  -Ij-(C„,  -1,-1.) 


x,x, 
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]j(c„  -xc„)dxdx  =  Lt(V2)C»,  -(I„  -I,)]-[c„,  -(I„  -I.)l 

X,  Xi 

])(<:(, 8e  -c„f)jxdx  =  L^(l/2)i.  -I,8j,  -  (I, 

XlX, 

~  [~ ^4  ~(^9  "^7)^1  “I4^E2  “(^10  “^8)^2] 
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Appendix  D.  Trim  Solution  Parameters 

bi=a, +ip(a2+a3)-p'(l,-lj 
Jet  =~(g|dXJ3T  V)J4t] 

JNi-Kl.-I.Xpd/V^^, 

JBF='g,(l.-I.)(p4/V)C™, 

Nt  =g2[llT  +-'2T(>P<i/V)]  +  p’S„ 

JjT  for  j  =  1,23/4,5,6  is  Jj  with  replaced  by  Sex.Tt 

E,E=E,+g2(xd/V>:E„ 

E3  =(Ei+ipE2L-p^x)/(0^ 

E4  =g4[cfi  +(ipd/V)c„  -Cd] 

E5  =-pVoo  +g4[cfiM/V)+^ 

dx  ^ 

ScT  =(VL)f8Erdx 

Nt  =  [1  +  ix(pd/V)]NT(iOo^ 

Egp  =[Egp  +  ix(pd/V)CpjEp]coQ 

Eb  'g.k,r,  +(c„r,  -c„8,,Xip<i/v)-%^l+(p/<i).)’8, 

dx 

fn  =I2(4t  -rT(x,))+[UT  +  W(x  J+(i4  -i,)rT(x,)Xipd/v) 

“Pt(xi))+[Ii34t  ■*‘^3^et(xi)+(I9  “ I? )Et (xj )Xip/^)“ ™idt 
nijT  =I4(^T  ~Et(Xj))+ [l,4^T  ■^I4^Et(x2)‘*'(IiO  ~  ^8  Xt  (Xl  )Xipd/V)  — Hljuj 

®^1DT  ~  IlD[^ET(Xl)~8Tc]+l3ijrj(x,) 

®^2DT  ~I2d[^Et(x2)“^Tc]‘*‘^4dEt(x2) 
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Appendix  E.  Transient  Solution  Parameters 


bn  =  a,  +ipa,  +  (a,  -ipI,)Aj  +  I.Aj 

=(AJ,  -ipI.X<i/vK 

Nw=gJJn+J.»(A^d/V)]-Ajvn 

dw. 

for  7  =  1,2,3, 4, 5,6  is  Jy  with  replaced  by  - - 

Ejk  =(E,+AkE2L+A^x)cOo" 

Edk  =g4[cfi  +(Akd/vK2  -Cd] 

Esk  =  Ak/(0^  +g4|^Cf,(Akd/V)+^j 
Efik  =  l  +  2k©;'’(Ak-ip) 

Vkc  =(VL)/v|/kdx 
N:=[l  +  x(A>d/V))Nj< 

4  =  i{*  + [i„  +  i,¥,(x, )+ (I,  -I,  )^;!;^j(A,d/v) 

+I=¥,(x,)+(I. -ijd¥^^^)j(A„d/V) 
m.k  =  ^l3N^k(x.)-.(l,  -I,)^](A,d/v)-m 

m2k  =  +  Ii4  +l4M'k(x2)+(lio  - Ig ) j(Atd/V ) - 1 

m,Dk  =  IlD[Vk(Xl)-M^kc]+l3D 
mjDk  =l2D[¥k(x2)-M^kc]+l4D 


f2k=l2  1 


m,k  =13  1- 
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Appendix  F.  Missile  Parameters 


Flare-Stabilized  Missile  Parameters 

L  =  25 
d  =  0.20  ft 
m  =  2.40  slug 

=0.0132slug-ft^  (0) 

I,  =5.006slug-ft^  (5.000) 


0.80ft 
S  =  (V4Xd)^ 

V  =  18000ft/sec 
p  =  .002slug/ft' 


o 

tl 

o 

-5<x^l2.5 

=  -1.14(x  +  5) 

-12.5<x<-5 

0 

K> 

II 

0 

Cn„=32 

Cm„  =-320 

Finned  Missile  Parameters 

L=20 

V  =  6000  ft/s 

d  =  0.35  ft. 

p  -  0.002  slugs 

m  =  3.50  slug 

^01  =^23  =1 

/,=  0.054  slug-ft2  (0) 

a,  =  1100  ft/s 

I,  =  14.318  slug-ft2  (14.292) 

1 

II 

c 

0 

10<x^ll 

_  4g7(x-10) 

-5<x:Sl0 

=  -(2/jtXl5  +  3x)aL 

-7<x^-5 

II 

r 

-ll<x<-7 

c„=2(ll-xX 

10<x^ll 

=  2 +  0.571(1 

-5<x^l0 

=  2.571+ (l/3jcXl5  +  3xXaL 

-7<x<-5 

=  2.571 -(l2/7tX6  +  x)aL 

-ll<x<-7 

Cd  =  (0.30X1 1  -  x)  -  0.295e'‘®<*‘‘®) 

10<x^ll 

=  0.0050 

-5<x^l0 

=  -0.10(5  +  x)+ 0.005  -  0.198e“*<’‘*’^ 

-7<x<-5 

=  0.0070 

-lUx^-7 

=0.14 

Qp=-18 
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1 

o 

11 

CM 

^22  =  -0.25x10"' 

CNa=9.7 

Cmc  =-34.4 

Cm.  =-980 

Cm.  =-190 

List  of  Symbols 


A(x) 

Cfj(x) 

d 

E 

A 

F. 

A 

Fd 

A 

Fe 


gl 

g2 

g3 

Sa 

S. 

St 

I 

Ix 

It 

A 

k 

L 

m 

Me 

P 

Q 

s 


local  cross-sectional  area 
aerodynamic  force  distributions  functions 

rod  diameter 
Young's  modulus 

complex  aerodynamic  shear  force  on  rod  at  x 

complex  beam  damping  shear  force  on  rod  at  x 

complex  beam  elastic  shear  force  on  beam  at  x 
pV^S/2 

gi/md 

pVSd/2m 

gaL/toJ 

Ix/md' 

I,/md^ 

(d)”*  JJy^dydz  =  (d)^  JJz^dydz,  area  moment  of  rod 

axial  moment  of  inertia  of  projectile 
transverse  moment  of  inertia  of  projectile 

beeim-damping  coefficient 
dimensionless  length  of  rod 
projectile  mass 

complex  elastic  moment  on  rod  at  x 
projectile  spin 

q  +  ir ,  complex  transverse  angular  velocity  of  body  1 
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V  magnitude  of  projectile  velocity 

location  of  beam  ends 


Xo,,X 


23 


P 

r 

6 

Sc 

Se 


dimensionless  length  of  fore  and  aft  aerod5mamic 
extensions 

angle  of  attack  of  central  disk 
angle  of  sideslip  of  central  disk 

— ,  cant  of  disk 

dx 

Sg-Sj  lateral  displacement  of  disk  relative  to  missile 
center  of  mass  (cm) 

L“*  [SecIx  /  lateral  location  of  missile's  cm 
*1 

bgy  +  ,  lateral  displacement  of  disk  relative  to  central 

disk 


roU  angle 


k 

K 

p 

Pi 

<j 

®0 

(Dl 

4 

F  =  (Fx,Fy,Fj 
H  =  (H„Hy,H,) 


frequency  of  k-th  mode 

damping  of  k-th  mode 
complex  cmgle  of  attack  of  disk 
air  density 

m/Ld ,  linear  density  of  rod 

®l/®R 

El/p,d^ 

lowest  elastic  frequency  of  beam  in  vacuum 
rigid  projectile  zero-spin  frequency 

+  la ,  complex  angle  of  attack  of  projectile 
aerodynamic  force  exerted  on  missile 
angular  momentum  of  missile 


50 


aerodynamic  moment  exerted  on  missile 
^  “  ('^x  >  )  missile  velocity 

^dc  ~  ('^xj'^dcy* '^dcz  )  velocity  of  disk  relative  to  missile's  cm 

dimensionless  location  of  disk  relative  to  central  disk 

dimensionless  location  of  disk  relative  to  missile's  cm 
angular  velocity  of  central  disk 
angular  velocity  of  non-spinning  elastic  coordinates 
reed  part  of  z 
imaginzuy  part  of  z 
Circumflex  superscript  denotes  a  beam  shear  force 
Tilde  superscript  denotes  elastic  parameter  for  bent  missile 
B  subscript  denotes  parameter  for  bent  projectile 
BF  subscript  denotes  bent  fin  parameter 
E  subscript  denotes  elastic  coordinates  at  central  disk 
T  subscript  denotes  trim  motion  parameter 


~  (^>^Ey  ’^Ez) 

S  =(x,8,,5.) 
a  =  (p,q,r) 

=  (0,q,r) 

•f{x} 


* 
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